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ABSTRACT 

Beam instabilities in the strongly magnetized electron-positron plasma of pulsar 
magnetospheres are considered. We analyze the resonance conditions and estimate the growth 
rates of the Cherenkov and cyclotron instabilities of the ordinary (O), extraordinary (X) and 
Alfven modes in two limiting regimes: kinetic and hydrodynamic. As a preliminary step, 
we reconsider wave dispersion and polarization properties in a one dimensional pair plasma 
taking into account relativistic thermal effects. We then find the location of the Cherenkov 
and cyclotron resonance of the X, O and Alfven waves with the particles from the primary 
beam. The importance of the different instabilities as a source of the coherent pulsar radiation 
generation is then estimated taking into account the angular dependence of the growth rates and 
the limitations on the length of the coherent wave-particle interaction imposed by the curvature 
of the magnetic field lines. We conclude, that in the pulsar magnetosphere the Cherenkov-type 
instabilities occur in the hydrodynamic regimes, while the cyclotron-type instabilities occur in 
the kinetic regime. We argue, that electromagnetic cyclotron-type instabilities on the X, O and 
probably Alfven waves are more likely to develop in the pulsar magnetosphere. 



1. Introduction 

At the moment the most promising theories of the pulsar radio emission generation are based on the 
plasma emission model, in which the high brightness radio emission is generated by plasma instabilities 



developing in the outflowing plasma (Melrose 1995). To find the instability that can be responsible for the 
generation of pulsar radio emission it is essential to know the dispersion relations of the normal modes of 
the medium and take into account the evolution of the modes as they propagate outward in the pulsar 
magnetosphere. 



In a standard model of pulsar magnetospheres (Goldreich & Julian 1969, 3turrock 196C) rotating, 
strongly magnetized neutron stars induce strong electric fields that pull the charges from their surfaces. 
Inside the closed field lines of the neutron star magnetosphere, a steady charge distribution established, 
compensating the induced electric field. On the open field lines, the neutron star generates a dense flow 
of relativistic electron-positron pairs penetrated by a highly relativistic electron or positron beam. This 
relativistic flow generates the observed pulsar radio emission. 

First, we present a consideration of the properties of linear waves in a strongly magnetized electron- 
positron plasma with similar distributions of electrons and positron. We take into account possible 
relativistic temperatures of plasma component. The particles are assumed to be in their ground gyrational 
state so that the plasma is one-dimensional. Properties of a one-dimensional pair plasma are considerably 



different from the properties of a well studied electron-ion plasma. Previously, Suvorov & Chugunov 1975 



investigated the dispersion relations in a one-dimensional plasma with a power law distribution of the 
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particles, Arons & Barnard 1986 considered waves in a cold plasma and in a warm plasma using a water 



bag distribution. For the propagation of the coupled O-Alfven modes, Arons Sz Barnard 198(: considered 



the case of infinitely strong magnetic field and found that the O branch is always superluminous and 



thus cannot be excited by the Cherenkov mechanism. On the other hand, Volokitin, Krasnoscl'skikh fc 



Machabeli 198E took into account finite magnetic field and found that the O wave becomes subluminous 



for large wave numbers allowing Cherenkov excitation. Lyutikov 1998 considered the waves in pair plasma 
taking into account the relative motion of electrons and positrons. Zank & Greaves 1995 used a fluid theory 
to consider waves in a nonrelativistic three-dimensional pair plasma. In our approach we use relativistic 
kinetic theory assuming that the possible resonant and nonresonant contributions from the beam particles 
may be considered as perturbations to the initial plasma state. 

We resort to similar distributions of the pair plasma with equal densities and neglect the curvature 
of the magnetic field lines. These approximations require some justifications. First, the pulsar plasma is 
nonneutral due to the presence of the primary beam. We neglect the nonneutrality of plasma since the 
density of plasma n p is thought to be much larger than the Goldreich- Julian density ugj — l~i • B/(2 7rec), 
rip = Xuqj = 10 3 — 10 6 n<3j (f2 is a rotational frequency of the neutron star, q is a charge of a particle, 
B is magnetic field, c is a speed of light, A is the multiplicity factor). Secondly, the inhomogeneity of the 
magnetic field results in a curvature drift of the particles perpendicular to the osculating plane. It can be 
shown that for the typical parameters in the pulsar magnetosphere the drift velocity of plasma particles 
due to the curvature of magnetic field lines could be neglected in the calculations of the real part of the 
dielectric tensor unless the curvature of field lines Rb satisfies the condition Rb <C 7p?\L (here 7 P is the 
average streaming energy of plasma particles in the pulsar frame and = c/lub ia a Larmor radius). This 
follows from the assumption of a nonrelativistic transverse motion and the transformation of the radius 
of curvature seen in the center of gyration frame Rb, eg = Rb/i 2 (Rb is the radius of curvature in the 
pulsar frame, 7 is the Lorentz factor of the particle). This condition is well satisfied inside the pulsar 
magnetosphere for the plasma streaming energy j p < 10 . By contrast, the drift velocity may be very 
important for the high energy resonant particles. 

In this paper we consider wave excitation in a strongly magnetized pair plasma in the approximation 



of straight magnetic field lines, thus omitting an important Cherenkov-drift resonance (Lyutikov, Blandford 



fc Machabeli 1998b| , [Lyutikov, Machabeli fc Blandford 1998a| ). This is an important mechanism that may 



be responsible for the generation of the cone type emission in pulsars. The electromagnetic Cherenkov-drift 
instability occurs in the kinetic regime on the high frequency vacuum-like O and X waves. It has the same 
advantages as the electromagnetic cyclotron instabilities considered in this paper. 

In presenting our results we were trying to find a balance between providing an exact general relation 
and finding a simple useful analytical approximation, which gives an idea of how modes behave in different 
regimes. Several approximations will be extensively used. For the case of radio waves propagating in the 
pulsar magnetosphere there is naturally a small parameter, u/lob (u> is a frequency of a wave), so that in 
many cases complicated dispersion relations and polarization properties can be simplified by expanding 
in this parameter. When making such expansion one should be especially careful near the points where 
dispersion curves almost intersect (Eq. |i"5j ). A formal expansion in lo/ujb diverges near the intersection 
point. In this case, since the intersection formally occurs only for the parallel propagation, the relevant 
small parameter is the angle of propagation with respect to magnetic field. Another simplification that 
we will often use, is the expansion in a ratio of the plasma frequency oj p to the cyclotron frequency 
ll>b = eB/mc (m is a mass of a particle) w p < wj. This is a good approximation for the most parts of the 
pulsar magnetosphere. 
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The main conlusions of our work are the following. For the chosen parameters of the magnetosphere 
plasma, the Cherenkov-type electrostatic beam instabilities develop in a hydrodynamic regime, while 
cyclotron-type electromagnetic instabilities develop in a kinetic regime. Electrostatic beam instabilities 
in the pulsar plasma are generally weaker than the electromagnetic instabilities. In addition, Chcrcnkov 
instabilities have largest growth rate near the stellar surface, where the Cherenkov resonance can occurs only 
on the Alfven mode. However, this mode cannot escape to infinity, even though it has some electromagnetic 
component. Another factor that limits the development of the Cherenkov-type instabilities is that they 
grow within a much narrower angles than cyclotron instabilities. In a curved magnetic field this results in a 
shorter length of the coherent wave-particle interaction. 

The relative weakness of electrostatic instabilities as compared to electromagnetic instabilities is 
an unusual characteristics of the strongly relativistic beams. The reason is that for the particles in 
the primary beam, which contribute to the development of the instability, the effective parallel mass is 
m e ffii = 7 3 m « 10 21 m. This suppresses the development of the electrostatic instabilities. In contrast, the 
effective transverse mass , m c g± = 71, m, is less affected by the large parallel momentum. The electromagnetic 
instabilities are less suppressed by the large streaming momenta. Thus, the relativistic velocities and 
one-dimensionality of the distribution function result in a strong suppression of the electrostatic instabilities 
as compared to electromagnetic instabilities. 

The calculations presented here provided a basis for the model of pulsar radio emission presented in 
Lyutikov, Blandford fc Machabeli 1998b. 



2. Plasma Parameters 

To a large extent a possible mechanism for the generation of pulsar radio emission is predicated on the 
choice of parameters of the plasma flow that is generated by a rotating neutron star. At this point we know 



only the general features of the distribution function of the particles in a pulsar magnetosphere ( Tademaru 



1973,Arons 1981, Daugherty fc Harding 1983] ). It is believed to comprise (see Fig. |l|) (i) a highly relativistic 



primary beam with the Lorentz factor 7^ rs 10 7 and density equal to the Goldreich- Julian density n G j, (ii) 
a secondary electron positron plasma with a bulk streaming Lorentz factor j p « 10— 1000, a similar scatter 
in energy T p « 7 P and a density much larger than the beam density n p w Xuqj = 10 3 — 10 6 ncj, (iii) a tail 
of plasma distribution with the energy up to 7 t = 10 4 — 10 5 . 

We will normalize the density of the pair plasma to the Goldreich-Julian density. 

n a = Xn GJ = 10 3 - 10 6 n GJ , u 2 p = \u 2 b = 2\v B tt (1) 

(subscript a in Eq. ([!]) refers to the electrons and positrons of the bulk plasma). Secondary pairs are 



born with almost the same energy in the avalanche-like process above in the polar cap (Arons 1983). The 
combination of the pair plasma and primary beam is expected to screen the rotationally induced electric 
field so that the flow is force-free. 

Another relation between the parameters of the plasma and the beam comes from the energy argument 
that the primary particles stop producing the pairs when the energy in the pair plasma becomes equal to 
the energy in the primary beam: 

2 < 7 >± n± — "f' b nGj, at the pair formation front (2) 



where < 7 >± and n± are the initial average energies and densities of pair plasma, 7^ is the energy of 
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the beam in the pulsar frame and uqj is the density of the beam. It is assumed that the initial densities, 
temperatures and velocities of the plasma components are equal. For cold components < 7 >± = j p , {p/ p 
is the stream 7-factor of the bulk plasma with respect to the pulsar frame), while for the relativistically 
hot components with a temperature T p the average energy is < 7 >jP — j p T p , where T p /2 is the average 
energy of particles in the plasma frame. 

In what follows, the quantities measured in the pulsar frame will be denoted with a prime. The 
relations between plasma parameters measured in the pulsar and plasma frames are 

7fc = 7^—1 Up = ——, uj b = uj b (3) 
L lp y/lp 

In this paper we neglect the difference of energies of secondary plasma components which arises as 
the flow propagates outward in curved magnetic field lines. For the consideration of the effects of relative 



velocity on the wave dispersion see Lyutikov 1998 



The uncertainty in the physics of the pair formation front forces us to allow for a broad range of plasma 
parameters. Accordingly, the growth rates of the particular instabilities can vary considerably depending 
on the assumed parameters. The numerical estimates will be given for a typical pulsar with the period P= 
0.5 s (light cylinder radius Ru — 2.4 x 10 9 cm), and the surface magnetic field B — 10 12 G and the primary 



beam Lorentz factor 76 = 2 x 10 7 (e.g., Arons 1983). These assumptions and the equation (g) reduce the 
number of free parameters to two: plasma temperature and the bulk streaming energy 7 P (or temperature 
and the multiplicity factor A) . Consequently, we will consider two separate cases of cold and relativistically 
hot plasma. For numerical estimates we will use the following fiducial numbers: 7 P = 100, A10 5 , T p <C 1 for 
the cold plasma, and 7 P = 100, A = 10 4 , T p as 10 for the relativistically hot plasma ( T p is the invariant 
temperature of plasma in units of mc 2 ). 

The radial dependence of the parameters is assumed to follow the dipole geometry of the magnetic 

field: 



Wfl(r) = lo b (Rns) 



( Rns \ 

V v J 



u p (r) = u> p {R NS )(?^ . (4) 



3. Response Tensor for a One Dimensional Plasma In Staright Magnetic Field 



In the limits of applicability of our simplifying assumptions, the dielectric tensor is (Lyutikov 



Machabeli fc Blandford 1998a) 



1 ;/.'/ 



. \ - 2 f d Pz fa \- ^pa f dp z 

Lu ^ J 7 3 n o2 L w 2 J 7 /a 

,2 



dp z £ (k 2 x + k 2 )v 2 z 



7 3 J]o2 
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a J I 



Here 



fi^ = w - fc z u 2 i w B 7 1 , = cj - k z v z , (6) 

where / Q are one dimensional distribution functions of the components a, v z is a velocity along the local 
magnetic field, 7 is a Lorentz factor of a particle, k x , k z and k y are the corresponding components of the 
wave vector and magnetic field is directed along the z axis. 

For stationary and spatially uniform plasma we can use Fourier analysis which reduces the problem to 
the following system of equations for the perturbations in the electric field: 

Aa/s£/9(w,k)=0 (7) 

where e Q ^(aj,k) is the dielectric tensor of the medium and 

A Q/3 = k a kp - k 2 c 2 5 a p + u> 2 e a p(uj, k) (8) 

The normal modes satisfy a dispersion relation 

Det\A af3 \=0 (9) 
whose roots determine the time behavior of the perturbations. 



4. Waves in Cold Pair Plasma in Rest Frame 

4.1. Dielectric Tensor 

In this section we consider waves in a cold, strongly magnetized, electron-positron plasma in its rest 
frame. If the average velocities of the electrons and positrons of the secondary plasma are the same, then 
retaining only nonresonant terms we find from Eq. (^|) the dielectric tensor for cold pair plasma with 
coincident distribution functions: 

CfTrff! 1 ~\~ n n €7/?/ 



B ~ 

2w„ 2 



e zz = 1 - 

^xy €yx ^xz ^zx £yz ^zy (10) 

where uj 2 — Aim p e 2 jm is the plasma frequency, lob — \e\B/mc is the nonrelativistic positively defined 
cyclotron frequency. 
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4.2. Dispersion of the Normal Modes 

Equation (||) with the dielectric tensor ( |l0| ) factorizes giving the three wave branches: X and two 
coupled O and Alfven branches with the index of refraction n given by 

2uj 2 

n 2 = 1 - 2 p 2 X mode (11) 

LOB 

uj 2 — 2uj p 2 J (uj 2 —ujb 2 —2uj p 2 ^ 



uj'* 



-a 2 2 — n — 2 2~To 2 2 Ta Alfven and O mode (12) 

j — uj uj b — 2 uj tjp -\-2 to b tjp cos y ' 

where 8 is the angle of propagtion with respect to the magnetic field (see Figures || and ||). 

Equation (^lj) describes the transverse X wave with the electric vector perpendicular to the k-B 
plane and equation ( |l2| ) describes the coupled longitudinal-transverse wave which has two branches: O 
quasi-transverse wave with the electric vector in the k-B plane and quasi-longitudinal Alfven wave with the 
electric vector along B. 



4.3. Parallel Propagation 

The normal modes of the plasma for the parallel propagation are given by 

uji = \plujp (13) 

n ' =1 + -4% (14) 
uj b - UJ t 

(subscripts / and t refer to longitudinal and transverse polarizations of waves) . 

For exactly parallel propagation, the dispersion curves for the O mode and Alfven mode intersect at 

uj* = V2lj p k*c m V2uj p (l + lu 2 p /ujb 2 ) , for £f < 1 (15) 

This intersection occurs only for the parallel propagation, while for oblique propagation the dispersion 
curves for O, X and Alfven modes are well separated. It follows that equation ( |l3| ) describes the O mode 
for k < k* and the Alfven mode for k > k* , while equation ([TJ) describes the X mode for all frequencies 
and the Alfven mode for k < k* and the O mode for k > k* . 

As the wave propagates in a curved magnetic field of the pulsar magnetosphere, it's path on the 



CMA diagram (see Sections 4.6) depends on the branch that the wave belongs to. In the linear regime 
in wave amplitude and with the adiabatically changing parameters of a medium, a wave always stays on 
a given branch. For example, original electrostatic wave emitted along the magnetic field can acquire 
electromagnetic components as it propagates in the presence of a curved magnetic field. The linear 
evolution of the longitudinal plasma wave, emitted originally along the field line, will drastically depend on 
which branch (O or Alfven) the wave actually belongs to. The propagation of O and Alfven waves in the 
inhomogeneous plasma of a pulsar magnetosphere differs considerably ( |Arons fc Barnard 1986 , Barnard) 



fc Arons 1986 ). If the plasma wave is emitted with the wave vector (or frequency) above the intersection 
point k > k* , then the linear transformation of the longitudinal plasma wave will follow the Alfven branch, 
which cannot escape magnetosphere. On the other hand, if the plasma wave is emitted with the wave vector 
below the intersection point k <,k*, then the linear transformation of the longitudinal plasma wave will 
follow the O branch, which may escape from plasma. 
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4.4. Oblique Propagation 

In the pulsar magnetosphere the waves that may be important for the generation of the observed radio 
emission have frequencies much less than the gyrofrequency. In what follows, we will often use the low 
frequency approximation when all the relevant frequencies are much less than the gyrofrequency. In the 
cold plasma in its rest frame this implies: ui <C ujb- The solution of equation (|ll]) in the low frequency limit 
describes a subluminous transverse electromagnetic wave: 



2u. 



2 



•- - k 2 c 2 [ 1 - ) = k 2 v\, u<^u B , Xmode (16) 



where va is the Alfven velocity in a strongly magnetized plasma. 

Solutions of equation ( fl2| ) are more complicated. The simple form for the dispersion relation may 
be obtained near the cross-over point, where the dispersion relation of the O mode crosses the vacuum 
dispersion relation or in the asymptotic regimes far from the cross-over point. 

Solving ( |l2|) with the refractive index set to unity we find the cross-over point for the O wave. 

ujI = k 2 c 2 = 2uj p 2 + ujb 2 sin 2 9 (17) 

Near the cross-over point, the approximate dispersion relation for the O mode may be found using the 
relation 



UJ — UJQ = — 



where K(u>, k) = is the dispersion equation for a;(k). 
FromEq. 

we find 

uj = ko c + K,(k — ko)c (19) 

where 



1 duj 

c dk 



ujb 4 cos 2 9 sin 2 i 



k—ko 



4uj p 4 + 2 uj b 2 uj p 2 sin 2 9 + ujb 4 cos 2 9 sin 2 < 



4 uj p a + ujb 4 cos 2 9 sin 2 I 



(20) 

where we used the assumption ujb S> uj p . From ( |20| ) it follows, that the behavior of the dispersion relation 
of the O wave near the cross-over point shows a very sensitive dependence on the angle of propagation. 
There exist a critical angle 9 C = 2uj 2 /uj 2 b at which the dispersion relation changes: 

uj= V2uj p + ^(k- ko)c, k 2 ^2uj 2 /c 2 if 9 < 2uj 2 p /uj B 2 

" = fe - 4 !^t"; 4 ^ 4 . fc o 2 c 2 =2^+ ujb 2 sin 2 9 if 9 » 2uj 2 /uj b 2 (21) 

For angles smaller than 9 C we can generally use the approximation of parallel propagation when considering 
the dispersion relations of the waves, while for larger angles we must take into account the effects of oblique 
propagation. 

The other limits when the dispersion relations for the O and Alfven waves may be obtained in closed 
form are the asymptotic limits far from the cross-over point. The large and small wave vector asymptotic 
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solutions arc 



UJ={ 



ijj = < 



k 2 c 2 1 



2 u) v cos 



2 u) p 2 sin 2 i 



2w p 2 cos 2 6 ( 1 



(" 



2 oj„ sin 2 oj„ sin 



Fc 2 



fc 2 c 2 cos 2 e ( i 



fc 2 c 2 sin 2 



wave 



if kc 3> oj v 



Alfven wave 



2 w p 2 + k 2 c 2 (l - fc2 ^T 2g ) sin 2 9 wave 



if kc <C w„ 



Alfven wave 



(22) 



(23) 



4.5. Infinite Magnetic Field 

In the limit of infinitely strong magnetic field the dispersion relations for the O (plus sign) and Alfven 
modes (minus sign) are ( Arons fc Barnard 198G| ) 



- u P ± 



v / fc 4 c 4 + 4w p 4 -4/c 2 c 2 Wp 2 cos(2( 



(24) 



The short and long wave length asymptotics are then given by (22) and (|23| ) with the magnetic field set to 
infinity. 

An important point in considering the wave excitation in the superstrong magnetic field is that we 
cannot neglect the very large but finite magnetic field. In the approximation of the infinitely strong 
magnetic field the O mode is always superluminous and thus cannot be excited by the Cherenkov-type 
resonant wave-particle interaction. In this limit any instability would occur on the Alfven waves which are 
strongly damped as they propagate out in the pulsar magnetosphere. When the finite magnetic field is 
taken into account, O wave becomes subluminous for the small angles of propagation and can be resonantly 
excited by the Cherenkov, cyclotron or Cherenkov-drift interaction with the fast particles. 



4.6. CMA Diagram for Cold Pair Plasma 



CMA diagrams (e.g. Budden 1985) are useful tools in considering wave propagation. It is a plot of the 
refractive index versus some functions of wave, plasma and cyclotron frequencies. We chose the following 
coordinates for CMA diagram: 

-Ma) 2 

where Y and X are the standard quantities in the magnetosonic theory. With this choice of coordinates 
the lines of constant Z are the lines of constant density and are independent of wave frequency. The lines 
of constant W are the lines of constant wave frequency and are independent of the density. The regions on 
the CMA diagram are separated by the resonance, where n — > oo, and cutoffs, where n — ► 0. 

Using Eqs ( |ll| ) and (12) we find resonance 

W = 1 X modes (26) 
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W = i + Z ± J (I + Z) 2 - 2Zcos 2 9 Ok Alfven modes (27) 



2 1 " V V2 

and reflection points 

VF = 1 + 2Z X modes (28) 

wTltv & Alfv6n modes ( 29 ) 

W — 2Z 

For the X wave the curve n = 1 corresponds to Z — (vacuum case). For the O mode n = 1 at 
VF = 2Z + sin 2 (cross-over point) and Z — (vacuum case). Other useful relations for the resonances of 
the coupled O and Alfven modes are 



\ + Z+ J{\ + Z) -2Zcos 2 6> = 




h + z-J{\ + z) 2 -2z C o^e= { or7 ^ 2Q (so) 



The CMA diagrams are plotted in Figs. || and ||. 



4.7. Polarization of waves in cold plasma 



To find the polarizations of the waves we construct a matrix of cofactors of A ( Melrose 197E ) : 



Aq,^ 71 k a kp Ti £ryry ~\~ S a k^y k^jC^y^j k a k~y€ry{3 A^/CyCfj^) 

1 



77 e 7'? e ')7 



(31) 



Then the polarization vectors may be chosen as columns of \ a /3. 
For cold plasma, the elements of X a /3 are 



Ax.t — I — 1 



2<4 



-LJ 2 + LOB 2 



-1 



2lo 



£ + n 2 sin 2 i 



A** = n 2 -1 + n 2 + 



A TO = - -1 + n" 



2w 



CiJ 2 — Wfl 2 



2w 



cos 6* sin = A ZK 



2w 



~W 2 + L0 B 2 



2 n 2 w B 2 w 2 cos 2 9 

LO 2 (-L0 2 + LOB 2 ) 



A zz = 1 - n 2 



2l0 2 



2lo 2 



, -L0 2 + LOB 2 I 

We note that these relations are exact in frequency. 



-LO 2 + LOB 2 



- n 2 cos 2 i 



(32) 
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For the X mode using (|llj) for the refractive index in (32) we find the polarization vector for the X 
mode ex = (0, 1, 0). For the O and Alfven modes using (|l2|) for the refractive index in ( |32"|) we obtain the 
ratio of the electric field components in the wave: 



E x _ (w 2 - ujb 1 ) (w 2 - 2 w 2 ) cot ( 
E~ z ~ w 2 (u 2 -lo b 2 -2lo 2 ) 



(33) 



For the point far from the cross-over point we can use the approximation of a very strong magnetic 
field to find 

E x 



-1 



2^ 2 



cot ( 



(34) 



Using the relation (B3h we can estimate the polarization of the O wave at the cross-over point. We find 



that 



E x ^ 2 ujp 2 I 
E~ z ~ wb 2 



(35) 



For the angles smaller than uj 2 /uj b the O wave is quasi-longitudinal at the cross-over point and for larger 
angles it is quasi-transverse. 

Relations and (|35|) allow us to find the normalized polarization vectors: 



eo 



eo 



eA 



jcos6> ( 



1 - 



,0,1 



2 uj v cos 



sine\ + o(-, 



{^.0.-1} 

J 1 n wo 2 esc Q sec 6 \ 

I 1 ' ' 5P J 



2^7 



2 cj p cos 
P 



(0 

{l:0,4sf}+o(! 



sinfl, 0, cos 9 ( 1 



kc ^ oj p 
/cc <C uj p 



(l_^)} + (j) * 



/cc <C w B 



(36) 



(37) 



(38) 



which are accurate to O 



5. Waves in a Cold Pair Plasma in Pulsar Frame 

In the pulsar magnetosphere the plasma is moving along the field lines with a bulk Lorenz factor 
7 P ~ 10 — 1000. In this frame the waves propagating in different directions have different dispersion 
relations. So that the dispersion equation for the X mode becomes a fourth order equation for w(k). To 
simplify the consideration we will use the low frequency approximation from the very beginning, i.e. we 
expand all the relevant relations in 1/wg. In what follows, the quantities measured in the pulsar frame will 
be denoted with primes. 
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For the forward propagating waves, which in the plasma frame has 9 -C 1, we obtain 



I 2 
V 



fc/ c ( 1 _ ^ + 7p " P 2 / n 



„ J- ^ 

^p' 1 > 372^ 

7 P 



k'ccosQ' 1 



,'2 



'7p 1 



c 2 fc' 2 sin 2 e' 
4 7 P Up' 2 



v^Wp cos 0' 
Upfc + 372 



if fc'c » 7 P Wp 
if fc'c <C 7 P Wp 

if k'c <C 7 P w p 
if fc'c > 7 p w p 



The plasma frequencies in the two frames are related by ui' p = Wp/s/%- 

The cross-over point, where the O mode becomes luminal, is 

o, .12 

LUq h 47 w B » 

7p 



(39) 



(40) 



6. Waves in a Relativistically Hot Pair Plasma 

6.1. Effects of Thermal Motion on Wave Dispersion 

In this section we consider wave propagation in the relativistically hot, strongly magnetized electron- 
positron plasma. The thermal motion of plasma particles affects considerably the dispersion of the Alfven 
mode at frequencies lu > uj p and the dispersion of the O mode frequencies uj oj p . Another important 
quantitative modification is in the dispersion relation of the X mode. An important factor for the excitation 
of the X mode is the difference of its phase speed and the speed of light. This difference is roughly 
proportional to < I/7 3 > (Eq. |38|). It is decreased considerably by the bulk streaming of the plasma. In the 
relativistically hot streaming plasma there are more particles with low Lorentz factors, that contribute to 
< 1 /7 s > , than in the cold plasma streaming with with the same average velocity. So for a given streaming 
velocity relativistically hot plasma has larger < I/7 3 > and larger growth rate. 



6.2. Distribution functions 

To estimate the thermal effects on the dispersion of the plasma mode we use the two kinds of 
distribution functions: (i) waterbag distribution 

/fe)=(r' (41) 
0, otherwise 

(here pr ~ mcrfT is the scatter in moments) and (ii) relativistic Maxwellian distribution (see also Appendices 
[b| and ^ for the calculations of the relevant moments of the distribution) 

ft 

fiPz) = 2K Jp T) ex P {-PtPvU»} (42) 
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here (3t — 1/Tp, T p is the invariant temperature, p M is a four- momentum of the particle, is four velocity 
of the reference frame, K\ is a modified Bessel function. In most of the calculations to follow we will assume 
that the plasma is very hot: p^ jmc 3> 1 and T p — y/p^,/(mc) 2 + 1 3> 1. 

Both these distributions are "fast falling" at large moments. This is an important factor for the 
dispersion relation of plasma waves (see below) . The advantage of the water bag distribution is that the 
various moments of the distribution can be easily calculated. The relativistic Maxwellian distribution is 
explicitly Lorentz- invariant (see Appendix [§| for details of Lorentz transformation) . 

The relevant moments of the distributions are summarized in Table [l] for the water bag distribution 
(in the plasma frame only) and in Table |^ for the relativistic Maxwellian distribution (in both plasma and 
pulsar frame). 

The water bag distribution is generally a good approximation for the account of the thermal motion of 
the particles. Its major drawback is the absence of a tail of high energy particles, that can resonate with 
the waves in the plasma. The Cherenkov resonance on the tail particles will result in a strong damping of 
the waves. The cyclotron resonance on the tail particles may result in a wave excitation if the distribution 
function is asymmetric with a long high energy tail. The condition, that the Cherenkov resonance is 
unimportant, is that the phase speed of the waves in plasma is much larger, than the thermal velocity of 
the particles. In the case of the idealized water bag distribution this condition has to be put in by hand. 
Whenever the phase speed of the wave becomes comparable to the thermal velocity the waves should be 
considered strongly damped and nonexistent. Therefore, we expect that the high frequency branch of the 
Alfven wave, which in the limit of cold plasma had a very low phase velocity, will be strongly damped. 

Here, we should also mention a long standing controversy about the dispersion of the longitudinal 
waves and the possibility of the two stream instabilities in the relativistic plasma. In the initial work ( 
Silin 1960| ) and later works ( Suvorov & Chugunov 1975) it was stated that the relativistic plasma does 



not support subluminous longitudinal waves. This problem has been considered anew (Tsytovich & Kaplan 



1972 ) who found subluminous waves. The controversy has been resolved by Lominadze fc Mikhailovskii 
197S who demonstrated the existence of the subluminous waves in the range 0<ti-l<l/<7> 2 (n~l), 
provided that the third moment of the distribution (< 7 3 >) is finite (here n is the refractive index and 
< 7 >S> 1 is the average Lorentz factor of the plasma particles). Thus, when the distribution function falls 
off at large momenta slower than \ subluminous plasma waves do not exist. 

For the water bag distribution, the dispersion of the plasma waves for the parallel propagation is given 
by Eq. (Eq). We find that n — 1 becomes larger than 1/ Tp for u > several times ujq. For larger frequencies 



the Longitudinal plasma waves are either strongly damped or do not exist at all ( Silin 1960) 



< 7 > 


< pv > 


<!> 

7 


< ^> 

7^ 


IT /2 


7t/2 


7T 


J_ 

7T 



Table 1: Relevant moments of the water bag distribution in its rest frame (dimensionless units). It is 
assumed that pT/m e c sa 3> 1. 





Plasma frame 


Pulsar frame 


Density, < 1 > 


1 


J_ 


< 7 > 


K Q +K ? , f (l + T„/2) if T„ < 1 
2 *"i 1 T p if T„ > 1 




7p (l + T p /2 + T p v 2 p ) if T p < 1 
7 P T P (l + «p 2 ) if T p > 1 


< v > 







<p> 





7p w p"fff — ' 


7 P ^(1+ T„/2) if Tp <C 1 
2 7pW pTp if T„ > 1 


< pv > 


T 


IP \ 2Ki V P ^ 1 P) 


7p(l+T p + ^) ifT„«l 
7pT p (l + t !p 2 ) ifT p »l 






(1 - T p /2) if T„ < 1 

V ifT ? » 1 


^ J 


ifT p « 1 
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6.3. Dispersion Relations in Relativistic Pair Plasma 

To simplify the analysis we will use the low frequency approximation lo <C lob and the assumption of a 
very strong magnetic field <C 1 from the very beginning. The dielectric tensor is then given by 

e xx = 1 + d T p (l + n 2 /3f, cos 2 6) = e yy 
2 n 2 lo 1 

£ " =1 T p (l-n 2 2 cos 2 9) +dTpn2 Sin2 9 

^xy = €yx — ^xz — ^zx — ^yz — ^zy (^^) 

where 

d= % h= f-^ (44) 

The normal modes of a hot plasma are given by the solution of (^) with the dielectric tensor (f43|). 
Similarly to the cold case, equation (^3|) factorizes into a dispersion relation for the X mode and a coupled 
equation for the Alfve and O modes. 



6.4. Parallel Propagation 

In the case of parallel propagation the dispersion equation gives two transverse wave with the dispersion 

to 2 = k 2 c 2 (1 - d T p (l + &)) (45) 

and a plasma wave 

u, 2 = ^ + fc 2 c 2 /3 2 (46) 

It is also useful to represent the dispersion relations for the plasma waves near the cross-over point in 
the form ([L9h. For the relativistic plasma components we find 



f dK{u,k) \ I ( dK{u,k) \ = f <j(l + v) 2 > 
uJ*kc-±(k-k ) (47) 



The phase speed for the high frequency asymptotic of the plasma wave (46) approaches the phase 
speed of the thermal particles c(3t- For the more realistic distribution function, these parts of the dispersion 
relation will be strongly damped on the Cherenkov resonance with the thermal tail particles. The high 
frequency asymptotic of the plasma wave belongs to the Alfvcn wave. From this we make a conclusion that 
the high frequency (kc > T p uj p ) part of the Alfven wave is strongly damped and does not propagate. 



6.5. Oblique Propagation 

The dispersion relation for the X mode is 

u? = k 2 c 2 (1 - dT p (1 + Pi cos 2 &)) (48) 
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The dispersion relations for the O and Alfven waves in a hot pair plasma are plotted in Fig. ^|. 

It is possible to obtain the asymptotic expansion of the dispersion relation of Alfven and O modes in 
the limits of very small and very large wave vectors. In the limit kc U T p ui p we have 



or = < 



■ 2 k 2 (1 - d T p (1 + pi cos 2 9)) ( 1 - 
while in the opposite limit kc <C T p lj p 



2 ui* sin 2 ( 



c 2 T p k 2 (-1+Pt cos 2 ( 



c 2 k 2 cos 2 9 (1 - d T p (1 + (3 2 cos 2 0)) ( 
^ + c 2 fc 2 (/3 2 cos 2 + sin 2 9) 



1 - 



c 2 fc 2 sin 2 I 
2 T„ w 2 



Alfven wave 
O-wave 



Alfven wave 
O-wave 



(49) 



(50) 



The X mode is always superluminous and the Alfven mode is always subluminous. The O mode is 
superluminous for small small vectors kc -C yfT v \o v and may become subluminous for very small angles of 
propagation 9 <C yjT p ~w p / Vb- 



The cross-over point (where the phase speed of the O mode become equal to the speed of light) is now 
= k 2 c 2 rs 2 TpUj 2 + lob 2 sin 2 9. Using relation (|l8| ) we can approximate the dispersion relation near the 



cross-over point as 



uj = kc — n{k — fco), 



1 (i4+47>f^ 2 )g 2 

y2 



16T 6 w 4 



(51) 



6.6. Polarization of waves in a hot plasma 



In the case of a hot plasma the matrix of cofactors A^g is quite complicated and is not given here. 
Simple relations may be obtained in the limit lob — oo and near the cross-over point for the O wave. In the 



limit lob = oo wc find the elements of the matrix A 



A« = (-Hn 2 ) -1- 



(h) 

a/3 



2 co 2 



T p uj 2 (~l + n 2 v 2 cos 2 9) 
\£) = n 2 (-1 + n 2 ) cos 9 sin (9 = \W 



2 • 2 

n sin t 



iW = 1 - n 2 - 



uu 



2lj 2 (-1 + n 2 cos 2 9) 
TpLo 2 (-l + n 2 v 2 cos 2 t 



A^ = (-1 + n 2 ) (-1 + n 2 cos 2 9) 
For the X mode we find that the polarization vector is ex = (0, 1, 0), while for the O mode 



Ej_ 

E, 



n cost) smf 
— 1 + n 2 cos 2 



(52) 



(53) 
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These relations are valid for the points not close to the cross-over point of the O wave (near the cross-over 
point the approximation ujb = oo is not applicable). Near the cross-over point, n = 1, we find 

E x ujb 2 

nr. = (54) 

For oblique propagation the behavior of the O mode at the cross-over point changes at 

9 w 4T p 3 d (55) 

For smaller angles the O mode is quasiparallel at the cross-over point while, for large angles, it is 
quasitransverse. 

The polarization vectors for the O and Alfven modes are then given by 

{cos e (l - ^#^) , 0, - ((l + 2 ^lr 2g ) sin 9)}, kc » u, p 
{^,0,-1} 6<: ^ tU = u W (56) 

lo 2 tan#" 

f^4 



Ah) 



Ah) 



1.0. 2T lj 2 < fcc«o> p (57) 



6.7. Dispersion Relation for Hot Pair Plasma in Pulsar Frame 

The dispersion relations for the forward propagating modes in the pulsar frame in the limit ui' <C 



arc 



w'v = k'c (l - 



Y 

* W B 'P 



Coo= k'c(l-^r iffc'c » v ^7 P ^and0'<l/T p 

Co A = fc'ccosfl' (l - if kc « y/T plj » p (58) 

These relationships are valid for the frequencies satisfying the inequality 

lu < -f p uj B /T p (59) 

This is a condition that in the reference frame of the plasma the frequency of the waves is much smaller 
that the typical cyclotron frequency of the particles uob/T p . 



The cross-over point in this frame is 



O, J2rp 

^ 2 = ^+4 7 >|0' 2 (60) 



7. Hydrodynamic and Kinetic Instabilities 



The description of the beam-plasma instabilities is based on the scheme used to solve the general 
problem of linear oscillations in plasma. The initial equations are the linearized kinetic equations for the 
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particles in a self-consistent electromagnetic field and Maxwell's equations. When the unperturbed state of 
the beam and plasma is stationary and spatially uniform, we can use Eq. (|^) to find the normal modes of a 
medium. 

For the beam-plasma system the dielectric tensor e Q ^(o>, k) may be represented as a sum of contributions 
from plasma and beam. 

/. , i-\ r i *"" C -plasma . ^ 1TC bearn \ 

e Ql s{uj, k) = d a( 3 + —v al3 + —<r a ( 61 ) 

where a ^ sma a nd cr^g m are the conductivity tensors of plasma and beam. 

Sometimes it is possible to consider beam as a weak perturbation to the system. Then, in the zeroth 
approximation, the normal modes of the medium will be determined from (|) with 

^o/?" 1 to zero. This 

will produce a set of normal modes of the medium |oj(k)'}. 

If the plasma alone is stable, then the frequency of the normal modes will have a zero imaginary 
part. In the first approximation, dispersion relation (^) may be expanded taking into account a small 
contribution to the dielectric tensor from the beam. The frequency shift A(k)' of the normal mode w(k) is 
then determined from 

+ K b {uj,k) = (62) 

w(k) ! 

where K p (lu, k) and Kf,(ui, k) are the plasma and beam parts of Eq.(|^). For stable plasma without a beam, 
_Kp(w,k) and {w(k)'} are real. 

Two separate cases may be distinguished here depending on whether the complex part of the beam 
contribution to the dispersion relations ( |62| ) -Kb(w,k) is zero or nonzero. If Im.(Kb(u>, k) = then equation 
( |62| ) has real coefficients. The complex solutions of Eq.(^) (if any) are complex conjugates. Solutions 
with the positive complex part correspond to the growing waves. These are hydrodynamic instabilities. In 
hydrodynamic instabilities, all the particles of the beam resonate with the normal mode of the plasma. 
This requires that the growth rate of the instability be greater than the intrinsic bandwidth of the growing 
waves: 



A(k)' 



|k-<Jv|<Im(A(k)). (63) 

Here k is the resonant wave vector, Sv is the scatter in the velocity of the beam particles. This is 
satisfied for a very small scatter in the velocity of the beam particles, so that all the particles from the 
beam resonate with the beam. 

Alternatively, if the complex part of the the beam contribution to the dispersion relations ( |62| ) Kb(u), k) 
is nonzero, the frequency shift A(k) will always have a complex part. If the complex part of A(k) ; is 
larger than zero, then the corresponding normal mode w(k) will be growing at the expense of the beam 
energy, while for negative A(k) the mode will be damped on the resonant particles of the beam. This case 
corresponds to the kinetic instability. The requirement that the frequency shift A(k) z due to the complex 
part of -Kb(w,k) dominates over the shift due to the large real part of -K&(w,k) requires that the growth 



rate be much less than the the intrinsic bandwidth of the growing waves (reversed inequality (63)). This is 
satisfied for a very large scatter in the velocity of the beam particles, so that at any given moment only a 
small fraction of the beam particles is in resonance with the wave. 



Though the physical interpretations of the kinetic and hydrodynamic instabilities are quite different, 
they may be considered as two limiting cases of a general beam instability. For a relativistic beam traveling 
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along a magnetic field with average Lorentz factor 7;,, scatter in parallel Lorentz factors A7, and average 
pitch angle ip the condition of the hydrodynamic approximation (|6^) takes the form 

fc » c (t + + k ^ J + f ^ « r (64) 

where s is the harmonic number (s = for Cherenkov resonance, s ^ for cyclotron resonance) and T is a 
growth rate of an instability. For the kinetic instability, this inequality is reversed. 

From ([m]) it follows that there exist a critical pitch angle 

ipcrit = -\— (65) 



For ij) > ipcrit the scatter in pitch angles dominates over longitudinal velocity spread. For tp > ip cr it the 
average "longitudinal" mass of the beam particles decreases by the factor of (ip'Yb) 2 so that the instabilities 
whose growth rate is inversely proportional to the "longitudinal" mass of the particles (like Cherenkov 
instability of plasma waves) may be enhanced considerably. 

Relativistic particles propagating along the curved magnetic field of a pulsar magnetosphere initially are 
in the ground quantum state (zero pitch angle) . They can develop a finite pitch angle by (i) particle-particle 
collisions, (ii) interaction with the electromagnetic field (Compton scattering on the diffuse thermal photons 
or recoil due the emission of electromagnetic waves at anomalous cyclotron resonance), and (iii) when the 
adiabatic approximation for the propagation breaks down (when the Larmor radius becomes comparable 
with the size of the inhomogeneity) . The pitch angle of the particles is then determined by the balance of 
these forces on one hand and radiation damping at the normal synchrotron resonance and the force due to 
the conservation of adiabatic invariant on the other hand. 

In magnetosphere, the particle-particle collision time is very long compared with the dynamical time 
because of the relatively low density of particles, high speed and the one-dimensional character of the 
motion. We also assume that the Compton scattering on the diffuse thermal photons is unimportant and 
that the adiabatic approximation for the propagation of particles is satisfied. The transverse component of 
the force due to the radiation damping at normal synchrotron resonance dominates the transverse motion 
of the particles near the neutron star, making the pitch angles equal to zero. Then the pitch angles will 
remain zero throughout the region where the above conditions are satisfied. 

In what follows we assume that plasma is one-dimensional before the development of instabilities. The 
condition of hydrodynamic approximation ( |64|) is then 

A7 sw B A7 

k W c — + — — ^ T ( 66 ) 
In the kinetic regime this inequality is reversed. 

For an instability to be important as a possible source of coherent emission generation, its growth rate, 
evaluated in the pulsar frame, should be much larger than the pulsar rotation frequency Q. The growth 
rates in the pulsar and plasma frames are related by 

T' = — . (67) 

7p 
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So the requirement of a fast growth in the plasma frame is 

» 1- (68) 

Another, more stringent requirement on the growth rate comes from the angular dependence of a 
growth rate. The emitting plasma propagates in a curved magnetic field. If an instability has a considerable 
growth inside a characteristic angle 59' , then the growth length should be larger than 58 R c , where R c is the 
curvature of the magnetic fields. In the plasma frame this requirement is 

where we used 59' « 50 /i P . 



8. Cold Pair Plasma: Resonances 

In the cold plasma approach the resonant interaction between the fast particles and the plasma may 
be considered as the interaction of the waves in the plasma with the waves in the beam. The interaction 
is the strongest when the dispersion relations of the waves intersect. Consequently, we are looking for the 
possible resonances between the waves in the plasma (|l^) and the waves in the beam (see Fig. and|): 

uj = Vbk z (70) 

LOB 

uj = Vbkcos9± — (71) 

76 

As we will see in Section |l0|, the resonant interaction of the plasma waves with the Cherenkov waves 
in the beam J70| ) is described by the cubic equation for the frequency shift, which always has complex 
conjugate solutions. This implies that the Cherenkov resonant interaction of the waves in the beam and in 
the plasma is always unstable. 

In contrast, the frequency shift due to the cyclotron interaction of the waves in the beam and in the 
plasma (fnj) is described by a quadratic equation, which has two real solutions for the plus sign in (0) and 
two complex solutions for the minus sign in (|7l]). Thus, only the minus sign in (|7l| ) will contribute to the 
instability growth rate. The resonance (^Tj) with the minus sign is called anomalous Doppler resonance. 
This corresponding instability may be considered as the interaction of the negative energy wave in the beam 
with the positive energy wave in plasma. Due to the resonant coupling, the amplitudes of both waves grow 
exponentially. 

Now let us consider the condition for the resonances (|7p| ) and ( |7l| ) to occur. From the low frequency 
asymptotics of the Alfven waves (^) we infer that the possibility of the Cherenkov excitation of the Alfven 
waves depends on the parameter 

(72) 



2j b uj p 



LOB 

If /i < 1, then Alfven waves can be excited by Cherenkov resonance.n However, if fi > 1 then Alfven 



2 In the case of cold plasma this may be considered as a sufficient condition for the Cherenkov excitation of Alfven waves. 
In the case of hot plasma this is only a necessary condition (see below). 
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waves cannot be excited by Cherenkov resonance. Instead, resonance can occur for an O mode subject to 
the requirement of sufficiently small angles of propagation Fig. |§|. 

For the cold plasma in the region of open field lines we have 

- .->"< - < ra > 

So, at small radii (fi -C 1) it is the Alfven wave that is excited by the Cherenkov resonance, while for 
larger radii (/z > 1) it is the O-mode that can be excited by the Cherenkov resonance. In the outer parts of 
magnetosphere (r > IOORns) the parameter \x becomes much larger than unity: fi ^> 1. 

For the parallel propagation (and only in this case) the parts of the O and Alfven modes that have 
longitudinal polarization may be considered as forming a single plasma wave with a dispersion u> = \/2u> p . 
In this particular case, the excitation of either O or Alfven part of the longitudinal plasma mode is 
very similar. But as the waves propagate in the curved magnetic field lines, the parts of the plasma 
mode corresponding to the O or Alfven wave will evolve differently resulting in a different observational 
characteristics of the emergent radiation. 

In what follows we consider separately the two possible cases of Cherenkov resonances: /i > I and 
/i < 1. 

We also note that the X wave cannot be excited by the Cherenkov resonance. Though the formal 
intersection of the Cherenkov wave in the beam ([70]) with the dispersion relation of the X mode is possible 
for all frequencies if fi = 1, the transverse polarization of the X mode excludes a resonant interaction with 
particles streaming along the magnetic field. 

The cyclotron resonance on the X modes occurs at u> Tes <C u>b, provided that 

^ » 1 (74) 

UJB 

Using the fiducial plasma parameters of the cold plasma, we find 

^ = , 7l J2_ = H«_ = 1 . 3x W M 3 (75) 

LOB "fpUB 7p W B \RNS J 

which implies that the X mode can be excited by the cyclotron resonance only in the outer parts of 
magnetosphere for radii satisfying 

[rns) > [x<n 



( > I zidi I ^ 2 H r (7 (>) 



The location of the cyclotron resonance on the X mode is quite sensitive to the choice of the bulk 
streaming energy. Comparing the resonant frequency (Table |^) with the plasma frequency, we find 



^res.X.cycl _ ^g 3 ~fp /2 /Wb\ 3 / 2 



> I (77) 



u p 7 6 w 3 7bA 3 / 2 V 20 / 

which implies that the X mode is always excited with the frequencies much larger than the plasma frequency. 
Locations of the resonances in the cold plasma are given in Table |. 



Resonances 


Extraordinary wave 


Ordinary wave 


Alfven wave 


Cherenkov 




^■res £ — ' 


4^ A « 1/76 

// > 1, A: res > Wp/c 

flwfl^fi e > i/ 76 

ss « + sin 2 + ^(1 + 7^) 
// > 1 


\/2lo p co&6 if < l/7 6 
0J rea = < , \x > 1 = oo) 
V2^cot9 if0>l/ 7 6 


Cyclotron 


Wres = if " ° 2 ) 
{>■ > 1 


k ^ c = (i + l^(i + ^ 2 )) 


^res C ' 


\/2uj v cos 0, A: res c = — ^ — 5 u < 1 

2 1 2 
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9. Hydrodynamic Instabilities In Cold Plasma 

9.1. Dielectric Tensor for Cold Beam-Plasma System 

The dielectric tensor for beam of the density rib propagating with the velocity Vb along the magnetic 
field B through a plasma of the density n can be found from a general expression (Eq |) with zero drift 
velocity u a = and distribution function f a (p z ) = n p 5(p z ) + rib5(p z — Pb) (nb is a density of a beam and 
p z is a momentum of beam particles): 

i Up ujb io 



■cob 2 7f, w 2 a; 2 yy 



—i ujb 2 
jb 2 0J 2 uj 2 yx 
kujb 2 w Vb sin 9 

9 TT^ — € Z x 

76 or u* 
i kub 2 ub Vb sin 9 

7b 2 w 2 w 2 zy 



_ _ 2 u p 2 _ u b 2 _ k 2 u b 2 v b 2 sin 2 9 

9 ^ "9 9 ~9 \ / 

u z ^b u 7b ' W 



where u = u — kvb cos9, u 2 = (u — kvb cos9) 2 — ub 2 /% and 7^ = 1/y 1 — 

We will always assume that beam can be considered as a weak perturbation, so that we can employ the 
expansion procedure described in Section [jj. 



9.2. Parallel Propagation 

In this section we calculate the growth rates for the beam instabilities for the waves propagating along 
the magnetic field, that we will use later as guide lines for the general case of oblique propagation. 

For the propagation along the magnetic field the dispersion relation (^|) with a dielectric tensor |78| 
factorizes: 



2ul ul 



-l + -i + -jL = (79) 



-l + " 2 + 5 P 2 + 5 7T —, - =0 (80) 

- u B z 7ft ^ (± ws/76 + u) 

Equation (|7^ ) describes hydrodynamic excitation of longitudinal plasma waves. As discussed above, 
this may be a longitudinal part of Alfven or O mode depending on the parameters of plasma. 

Equation ( |80| ) describes the cyclotron excitation of the O and X modes. For the parallel propagation 
the cyclotron excitation of the Alfven wave does not occur. 
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9.3. Cherenkov Excitation of Plasma Waves for 6 = 



We now look for the correction to the relations (|T^) and (|7(]) when the two intersect. 



oj = 
uj = 



V2ujp + A 



Vb k cos 6 + A 

Expanding in small A, we find that the frequency shift satisfies a third order equation: 





(81) 



V2A 3 u b 2 



lb 



(82) 



Equation (|82j) always has one real and two complex conjugated roots. The complex root with the 
positive complex part corresponds to the instability. 



Solving Eq. (|82|), we find the complex part of the frequency shift: 



Im( A) 



2 e 7 & 



2 2/3 



7&V> 



(83) 



This is a growth rate for the Cherenkov excitation of plasma waves (c.f. podfrcy et al. 1975 , Egorcnkov 
et al. 1983|) . 



We can compare the importance of the Cherenkov excitation of plasma waves by evaluating growth 
rate ([83]) for the set of fiducial parameters of a cold plasma and comparing it with the dynamical time (p7u: 



Im( A) _ A 1 / 6 
7 P 



3/2 

Iblp 



86 



Rns 



-3/2 



(84) 



From which it follows that this instability may be important for r < 20. We will see in Section fL0| that 



the second criterion (69) is not satisfied for the Cherenkov excitation of Alfven or O waves, so that the 
Cherenkov instability does not develop. 



9.4. Cyclotron Excitation of Transverse Waves for 6 = 

We expect that the hydrodynamic instability will be strongest at small wave vectors. We can then use 
the low frequency approximation (|l6|) to the dispersion of the transverse waves. We seek the correction to 



the relations (16) and ( |7l| ) when the two intersect. 

,2 



kv b cos6 + A (85) 

7b 



Expanding in small A, we find that the frequency shift satisfies the quadratic equation: 

_A^ ± W^ Q 

Ld a A 7f, C<J Z 
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The ± sign in (|3q) corresponds to the two signs in (BOJ). For the normal Dopplcr resonance (plus sign 



in (80) and fl86h) the resulting frequency shift is real. For the anomalous Doppler resonance, the frequency 



shift is complex: 



= * \^ hkc ( 87 ) 



which gives near the resonant frequency (Table 

l Up LUf, 



Im(A) = 



For the parameters of a cold plasma the growth rate ( J8q ) is much longer than the dynamical time 
everywhere inside the light cylinder. 

Im( A) a/A a , , 

Which implies that hydrodynamic regime of the cyclotron instability is unimportant. 



9.5. Perpendicular Propagation 

Next, we consider the hydrodynamic instabilities for the waves propagating perpendicular to the 
magnetic field (magnetized Wiebel instability (Weibel 1959)). The normal modes of plasma without a beam 
for 9 = tt/2 follow from Eq. (g): 

LU 2 = k 2 +2LU p 2 (90) 

uj 2 =u B 2 + 2u p 2 (91) 
n 2 = 1 + 1^ - (92) 



In the limit ojb — > oo the two solution of the biquadratic Eq. fl92|) may be expanded in large ojb'- 



uj z =oj b z + 2oj p z + z-e- (93) 

OJB 



w 2 = c 2 fe 2 ] (!14) 



In the limit he <*C lob the dispersion curves ( |93"| ) and (^T|) approach each other near the upper hybrid 
frequency: lo 2 = lob 2 + 2oj p 2 . We then expand determinant J78| ) for 9 = tt/2 near the upper hybrid 



frequency oj — \/ojb 2 + 2oj 2 + A keeping the terms up to the second order in A. 



OJ^ UJ B to* 7& ub 

Solving Eq. (pa), we find the frequency shift 



A 2 ujb 2 , c 2 Ak 2 c 2 k 2 oj b 2 



kuj p 2 (Jjbk - \/j b k 2 - Auj h 2 ) 
A = " o r- 3 ( 96 ) 
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which shows an instability for kc < 2w(,/7& with a maximum growth rate 

. x k>r> LOh 

Im( A) mM » p 3 \ (97) 
which is negligible for all reasonable pulsar plasma parameters. 



10. Oblique Wave Excitation in Cold Plasma in the Hydrodynamic Regime 

In this section we develop a general theory of the hydrodynamic weak beam instabilities in the 
cold magnetized electron-positron plasma. We expand ( |78| ) in small u>\, keeping only first terms. After 
considerable algebra we obtain 



lo 2 + lob 2 J VV lo 2 J \ lo 2 —lob 2 

2 , , 2 , , 2 „„ 2 , 



2n u>b tjj p cos 

UJ 2 (LU 2 — LOB 2 ) 



3, .,2 \ w'-Wr' / \ L0 Z — Ct>R Z 



ry*LJ 2 \ IV -LOB J \ LO^-LOB 

LOb 



2 (2 ken 2 (-1 + n 2 + -4^-^) u) f3 b cos sin 2 



7b w 2 



/fc 2 c 2 (-1 + n 2 ) fj h 2 (-1 + n 2 cos 2 6>) sin 2 (9 



( 2{ -^ +UB2+ "fh n \?:^ ) sin 2 , 

P I W 2 (-CJ 2 +CJ B 2 ) 2 L0 2 {lO 2 -LOB 2 ) I 



to 2 I I 2ijJ, 2 \ / „ 4w c 2 , , 

p 1 ' -2 + n 2 + „ p „ + n 2 cos 2 i 
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The term containing 1/w 2 contribute to Cherenkov excitation and the term containing 1/lo 2 contribute 
to the cyclotron excitation. 



To find the growth rates we expand the plasma part of (£38|) near the plasma modes (Eqs. ([lT]) and 
(|l2|)) and the beam part near the resonances Co — (for Cherenkov excitation) or Co = (for cyclotron 
instability). The expansion of the plasma part of ( |98|) near the plasma modes is done according to the 
relation 



where K p is the plasma part of the determinant (|7 

2u. 2 



K-n = [1-n 



2 | " UJ P 



—UJ 2 + LOb' 
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LOB 2 J <-0 2 (w 2 -(i)fl 2 ) 



Simultaneously, in the beam part of the Eq. (|98|) we should use the normal modes of the medium for 
the estimates of u> and the refractive index n. 



and a/ ) are the solutions of the equation K p = 0. 



10.1. 



Excitation of Alfven waves 



We recall that the Cherenkov excitation of Alfven waves is possible only for [i < 1 (Section g). 
Expanding in small A near the Alfven wave dispersion relation in the limit of infinite magnetic field and 
using the resonant frequency (Table H) we can find the growth rate 



A = i 



~7 



2" 



7b 



V3 uj^ uli'S cot ( 



2E 7 , 



if « l/ 76 



if » l/ 76 



(101) 



The hydrodynamic growth rate of the Alfven waves has a maximum for parallel propagation. The 



numerical estimate of the maximum rate for Cherenkov excitation of Alfven wave is given by (84) subject 
to the condition \i < 1 ( Section |J) . 



We can now compare the instability growth length with the coherence length 
j I/76 we have 

R c S9lmA R c il A 1 / 6 [l7b~ 



Hp 



c 7p7 fc 2 V 

which implies that the Cherenkov instability on the Alfven waves is unimportant. 



< 1 



Using the estimate 
(102) 



10.1.1. Cyclotron excitation of Alfven waves 



Near the cyclotron resonance u) = we keep in ( Bq ) only the terms proportional to l/£j 2 . We consider 
two separate case depending on whether [i is larger or smaller than unity. 

For \i < 1 (i.e. near the neutron star surface) using the short wave length asymptotics of the Alfven 
branch (E2T) and the resonance condition (|7T[) we find a pair of complex solutions 



y/uj uj b V cos 6 sin A 1/4 lVcos6> sin6> /w B \ 1 /4 
A = ±l ; = TTU 1 77 J ( 103 ) 



This frequency shift agrees with the S_ of the Table 3, region 3 of Godfrey Godfrey et al. 1975. 
The ratio of the growth rate in equation ( |l03| ) to the dynamical time is 

Im(A) _ 76 1/4 \/Hoi0sin6» /w B y/4 f - \ ~ 3/4 



, < 150 - — (104) 
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The criteria of a short growth length (pq) gives 



R c 89Im(A) R C Q A 1 / 4 /h) fl y/* = / r \ 



-3/4 



(105) 



7p 

where we have approximated the trigonometric functions by unity. 



The equations (104) and (105) imply that the cyclotron excitation of Alfven wave in the hydrodynamic 
regime may be marginally important near the neutron star (see also Sections 13.1.2 and Eq. 127). 

For ^> 1 and 8 ^> (this case corresponds to the resonant frequency much smaller than oj p ) 

using the short wave length asymptotics of the Alfven branch ( p2] ) and the resonance condition ( |7l| ) we find 
a pair of complex solutions 



lu b ^ cj 6 tan(9 3 fitani /3 6> /cj b \ 2 / 3 

= ±l 2* 7 * W * ? 2W^p^ KIT) 



(106) 



Comparing this growth rate with the dynamical time we obtain 



Im(A) 



tan6»i/ 3 



( <£b 
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2/3 



3 x 10 4 tan6> 1/3 



(— 
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(107) 



where we estimated the typical angle of emission by unity. The ratio ( 107 ) is larger than unity for 
r/R NS > 43. 

The second criteria on the growth rate (|69|) becomes 



R c S6lm(A) R c n 



C7p 



1 /UJB_\ 2/3 _ 



c Ai/ 6 7 b 5/6 7p 7/3 v « 



= 120 



Rns 



(108) 



which implies that the cyclotron excitation of Alfven waves is efficient in the region where [i > 1 (^— > 43). 



10.2. Excitation of the X wave 

Using (^) we see that X mode is not emitted by the Cherenkov resonance since the corresponding 
term is zero if evaluated at the X mode (jl^) . 

Near the cyclotron resonance we obtain for the X mode a pair of complex solutions 

A = ± VBJ. 1Q9 

2 7b V w 



where we used kc w. Using Table ^ we obtain a growth rate for the cyclotron excitation of the X mode 



A = ±- — ^— — (110) 



76 w B 

The growth rate is almost constant inside the cone 8 < and is zero for larger 8. The growth rate 
(110) was also obtained for the case of parallel propagation (|8S|). As we saw in Section 9.2 the hydrodynamic 
instability on the X mode is unimportant. 
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10.3. Excitation of the O mode 



10.3.1. Cherenkov resonance ( (J, > 1) 



Starting with ( pq ) we expand the plasma part near the dispersion relation for O mode ( |l2| ) and the 
beam part near Cherenkov resonance ( |70| ) n = 1/ (/3& cos 0) . In the low frequency approximation uj <C u>b O 
wave can be excited by Cherenkov resonance only if /i > 1. We distinguish two cases: /J, > 1 and /i> 1. In 
the former case the resonance occurs at approximately cross-over point for the O wave while in the latter 
case the resonance occurs at u> ^S> loq. 



In the case fx > 1 with o> rcs given in Table |3[ the growth rate is 

y/E lob 3 oJpOJb 1 sin 9 



A = 



2 7b l (-1 + M 2 ) 7 



(8< 



2w B 2 sin 2 



(111) 



The growth rate is proportional to the angle of propagation with respect to the magnetic field. This is 
due to the increase of the potential part of the O mode with the angle. The case of fi — 1 corresponds to 



the unlikely case when the beam velocity is exactly equal to the Alfven velocity. The growth rate (111) is 
valid for 9 < I/7&. For all practical purposes this growth rate is very small. 



In the case /1 > 1 the complex shifts of the frequencies near the Cherenkov resonance to — A are given 



by 



Im(A) = 



V3t 



>p iOb' 
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lob 1 sin 2 1 
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2 7iUB 
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if 9 > 



2^ 

LO B 2 
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(112) 



The maximum growth rate is reached for the parallel propagation, while for oblique propagation the 
growth rate falls off as 9~ 2 I Z . We have already estimated the growth rate for the parallel propagation ( |84| ) 
and found that it can be marginally efficient. We can apply the second requirement on the growth rate 

which relates the characteristic angle of emission and the radius of curvature (J69|). Estimating 
find 

RcdeimA f\ i/3 \ / n\ fR c n\ 

(113) 



we 



C7p 



Jp 



< 1 



k lo b J \ c 

which implies that the Cherenkov instability on the O mode does not develop. 



10.3.2. Cyclotron excitation of the mode 

The growth rate for the cyclotron excitation of the O mode in the hydrodynamic regime may be 
estimated from the growth rate for the cyclotron excitation of the X mode (109) with the resonant frequency 
given m Table g. The maximum growth rate is reached for the parallel propagation and is equal to the 
cyclotron growth rate of the X mode (|109|). The growth rate (10E) decreases with the increase of the 
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resonant frequency. Since the cyclotron resonance on the O mode happens at the frequencies larger than 
the the cyclotron resonance on the X mode, the corresponding growth rate of the O mode are smaller for 
oblique propagation. 

We conclude this section by the table of the hydrodynamic growth rate in a cold plasma (Table Eh. 



11. Relativistic Pair Plasma: Resonances 

From the low frequency approximation to the Alfven waves dispersion we find that the possibility of 
Cherenkov excitation of Alfven wave in a relativistic hot plasma depends on the parameter 



(J-h = \ — 2 ~ 114) 

V u B 2 wb 



(compare with (]72 

Using our fiducial numbers parameter fi may be estimated 



ITXn ( r \ 3/2 <1, if hr- < 43 

^ = 2 7 Ji^^5xlCr 3 = \ R » S { (115) 

V ^Blp \RnsJ I > 1, if [^-j > 43 

Numerically /i^ and \i are equal for the chosen set of the fiducial numbers for the cold and hot cases. 

Similarly to the cold case, the parameter fi^ determines the possibility of the excitation of the Alfven 
and O waves. If /i^ < 1 then the O wave cannot be excited by Cherenkov resonance. In this case the Alfven 
wave may be excited by the Cherenkov interaction subject to the condition that the resonance occurs on 
the parts of the dispersion curve that are not strongly damped (see below). If /i/j > 1 then the O wave may 



excited by Cherenkov resonance for the angles of propagation 9 < 

Another limitation on the possible resonance comes from the requirement that the waves in the plasma 
are not strongly damped at the location of the resonance. This is an important constraint on the resonance 
of the Alfven wave, which is strongly damped at large wave vectors. 

Using the dispersion relation for the Alfven waves in the limit kc <ti lo p , we find that the cyclotron 
resonance on the Alfven wave occurs at kc <C lo p for the angles of propagation larger than 

lbto p 

For smaller angles the location of the cyclotron resonance on the Alfven wave depends on the parameter 

V = 3 (117) 

T p 2 uj b 

If r] <C 1 (very hot plasma), then the cyclotron resonance on the Alfven wave occurs in the region u lj^ \ 
where Alfven waves are strongly damped. If, on the other hand, rj 3> 1 (warm plasma), the cyclotron 
resonance on the Alfven wave occurs at approximately 0JQ h \ wher Alfven waves are not damped (Fig. ||). 

Since Alfven wave cannot escape to infinity, they should be converted to electromagnetic modes before 
they are damped on the thermal particles. The Alfven waves with large angles, which are generated with 
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the frequency uj -C uJq would have more time for the nonlinear processes to convert them into escaping 
modes, than the Alfven waves generated in a warm plasma with uj rs lu^ and small angles of propagation. 
Thus, the cyclotron resonance on the Alfven wave is likely to produce waves propagating in a cone around 
magnetic field. 

The resonances in the relativistic pair plasma are given in Table |[ 



12. Hydrodynamic Wave Excitation in Relativistic Pair Plasma 

12.1. Dielectric Tensor for the Beam-Hot Plasma System 

To simplify the analysis we will use the low frequency approximation uj -C ujb and the assumption of a 
very strong magnetic field -^r <C 1 from the very beginning. The dielectric tensor is then given by 

2-2 

( .,,,. 1 • d 7,. ( 1 • ,r ij cos 2 6) " '" '~ - 
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2n 2 u 2 p lrlTr 2-2n "I k 2 u b 2 p b 2 sin 2 u 
T p (l~n 2 (3 2 cos 2 0) +dTpn Sm 9 ~t£& jb^u; 2 (118) 



12.2. Parallel Propagation 

For parallel propagation Eq. with the dielectric tensor ( |11§| ) factorizes: 

2 n 2 uj? 

1 1— 5-r ^=0 (119) 

l-n 2 + „ P- % + d T p (1 + n 2 Pi) = (120) 

2 * 

l-n 2 „ f h f - + dT p (l + n 2 p 2 T ) =0 (121) 

7b U) Z (LOB/lb+UJ) 

Following the same procedure of expanding the dispersion relations in small frequency shifts A near 



the intersection of the two resonant curves, we find from ( 119 ) the growth rate for the Cherenkov excitation 
of plasma waves: 



2 6 lb V t p 2 ' lb y/lp T p 



(cf . with Egorenkov et al. ) . 
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Using the relations between parameters of the hot plasma (Eq. || with < 7 >= 2T P 7 P ), the condition 
of a fast growth (|6q ) for the growth rate (122) takes the form 



A V6 



Im(A) ^ 



3/2 ^1/2 v ^ I 



-3/2 



(123) 



For the fixed values of jb and 7 P the growth rate for the Cherenkov excitation of plasma waves in a hot 

2/3 

plasma is smaller by the factor T p as compared with the cold plasma. 



Solving Eq. (|12l|), we find the growth rate for the cyclotron excitation of transverse waves 



A = 



Tp uip uib 



'XTp n 



Comparison of this growth rate with the dynamical time gives 



(124) 
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(125) 



From (125) and (89) if follows that the cyclotron excitation of the transverse waves in the hydrodynamic 



regime is not affected by the relativistic temperature of the plasma particles and is not important in the 
pulsar magnetosphere. 

Similarly to the cold case we omit the details of the calculations of the growth rates and conclude this 
section by the table of the hydrodynamic growth rates in the relativistic hot pair plasma Table ^. 



13. Excitation of Oblique Waves in Relativistic Pair Plasma 
13.1. Excitation of Alfven Waves 

13.1.1. Cherenkov excitation of the Alfven mode ( /x/j < lj 

We give here the growth rates for the Cherenkov excitation of the Alfven mode in the limit /ij « 1. 
Then we can use the infinite magnetic field approximation to the dispersion of Alfven waves. We then 
expand the plasma part of the the dielectric tensor ( [43] ) in large gyrofrequency keeping the zeroth order ( 
d = 0). 

The complex part of the frequency shift is 

1 A 

\/3ujg u>b 3 cos 1/3 6 

Im(A) = 
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(126) 
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13.1.2. Cyclotron excitation of Alfven wave 



Following the same procedure as in Section 13.1.1 we find equation governing the cyclotron excitation 



of Alfven wave in the relativistic pair plasma in the case when the plasma part of the dispersion relation 
for the Alfven wave is calculated in the infinite magnetic field limit. Assuming that 1 — n cos 9 3> "f^ 2 we 
can set /3b to unity. 

Following the discussion in Section |ll| we are interested in the Alfven wave cyclotron instabilities 
occurring at frequencies u> < WTpU) p since for larger frequencies the Alfven wave is strongly damped. This is 
satisfied only for warm plasma ( r/ > 1). We can analytically find the growth rate in this case for the angles 
of propagation larger than /i^ (/i/,. <C 1). In this case, using uj rs fcccos#, we find two complex solutions for 
A 

, _ LUB^Ub tans 6 _ fitan 1 / 3 ^ /^b\ 2 / 3 , , 

Using the relations between parameters of the hot plasma we conclude that the cyclotron excitation of 
Alfven wave is unaffected by the relativistic thermal spread of the plasma particles. 



13.2. Cyclotron excitation of the X mode 

Following the same procedure and using the resonant condition (Table |5|), we find the following 
frequency shifts describing the cyclotron excitation of the X mode in a relativistic pair plasma: 

A = iV^v^ (128) 



13.3. Excitation of the O mode 

13.3.1. Cherenkov excitation of mode (/ih > lj 

Like in the case of the plasma, the Cherenkov excitation of the O mode case occur in two different 
regimes: [ih > 1 and /i/, 1. In the former case the resonant frequency is much larger that the cross-over 
frequency. In this limit the O mode is quasi transverse and Cherenkov excitation is unimportant. 

In the case [ih ^> 1 the Cherenkov resonance occurs approximately at the cross-over point. The 
resonant frequency may be approximated by the cross-over point: ui 2 = 2T p lo 2 + ui B sin 2 9. The complex 
part of the frequency shift is given by 
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(129) 



7b ljb 3 sin 6 



The maximum growth rate is reached for the parallel propagation (51). In Section 



10.3.1 



we found 
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that in cold plasma the Cherenkov growth of the O mode is unimportant due to the very short coherence 
length, which is, in turn, limited by the small range of angles of the growing waves. Since the hydrodynamic 
growth rate of the O mode in hot plasma is smaller than in cold plasma, we can make a conclusion that this 
instability is unimportant. 



13.3.2. Cyclotron excitation of the mode 

Similarly to the case of a cold plasma, the growth rate for the cyclotron excitation of the O mode may 
be estimated from the growth rate of the X mode with the resonant frequency given in Table H The growth 
rate has a maximum for the parallel propagation and decreases with the angle due to the sharp increase of 
the resonant frequency. 

We conclude this section by the table of the hydrodynamic growth rates in the relativistically hot pair 
plasma (Table ||) . 



14. Kinetic Instabilities 



As we have discussed in Section 0, a general beam instability may be treated analytically in the 
hydrodynamic and kinetic limiting cases. We have considered hydrodynamic beam instabilities in pair 



plasma in Sections 10 and |12|. Now we turn to the kinetic regime of instabilities. The condition for the 
kinetic consideration to apply is the opposite of the condition ( |63| ) . It requires a substantial scatter in the 
velocities of the resonant particles. In what follows we assume that distribution of the beam particles is 
described by the relativistic, one-dimensional Maxwellian distribution: 

1 I P,W \ 



^"'^t^^I-ttJ (130) 

where n\, is the density of the beam measured in the laboratory frame (the Lorentz invariant proper density 
is nt,jb), — (7b,/3b7b) is the four velocity of the rest frame of the beam, Ty, is the beam temperature in 
units mc 2 , K\ is a modified Bessel function. 

This function may be simplified in the limit of cold beam (in its frame) Tf, <C 1 and large streaming 
velocity 7& ^> 1. We find then 

/( ^7s exp l-^H (131) 

where pf = 7 2 T(,toc is the scatter in parallel moments. 



In case of kinetic instabilities the growth rate is given by (e.g., Melrose 197S ) 



(132) 

uj— aj(k) 



where e' aj3 and e'^ g are hermitian and antihermitian parts of the dielectric tensor, u>(k) is the frequency of 
the excited normal modes of the medium, and e a is its polarization vector. The antihermitian parts of the 
dielectric tensor are due to the resonant interaction of the particles from the beam at Cherenkov (|7(i| ) and 
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cyclotron resonances ([71]). Using the Plemelj formula we find 
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(133) 



Using the polarization vectors (|56|) , (j57[) we find that for the quasitransverse waves (O mode to ^> lo^\ 



Alfven mode lo <C lo^ and O mode lo ~ oOq 1 ' , 8 ^> uo 2 B / (T p lo 2 )) , while for the O mode at the cross-over 
point and 8 <§; lo b /(T p lo 2 ) 
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(134) 



With polarization vectors (|36|) and (|38|) we find from ( 133 ) , that for quasitransverse parts of the waves 
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(137) 



where we split the antihermitian part for the O and Alfven modes in two parts: e" ch is due to the 
Cherenkov resonance and e is due to the cyclotron resonance. 

Most of the relations (137), excepting e'^ ch , are valid for both cold and hot plasma. For hot plasma we 

have 
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(138) 



For the Cherenkov excitation of the O mode in the limit /i/j 3> 1 (when the resonance occurs at the 
cross-over point) we find 
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(139) 



The calculations of the integrals in (135 - 139) are given in Appendix [A| 
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14.1. Parallel Propagation 

We first consider an important, separate case of parallel propagation. 

Using the polarization vectors ei = (0,0, 1) for longitudinal waves and et = (1,0,0) for transverse 
waves we find 
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(140) 
(141) 
(142) 

(143) 



The corresponding growth rates are 



r* = 



Cf) r 



Tpkcw 2 \ d-f 



(144) 
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With the distribution function of the form ( |13l| ) we find 
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(146) 
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The kinetic growth rates ( J14q ) and (147) can be compared with growth rates in hydrodynamic 
regime (Eqns ( |85| ) and (]87|)). In a hydrodynamic regime both cyclotron and Cherenkov growth rates are 
proportional to the negative powers of the particle's Lorentz factor. This is a significant factor for the 
primary beam and for the particles from the tail of plasma distributions. In contrast, kinetic growth rates 



(146) and (147) are not suppressed by the relativistic streaming of resonant particles. On the other hand, 



kinetic growth rates (146) and ( |147j ) scale linearly with a small ratio of the beam density to plasma density 
while hydrodynamic growth rates (|83|) and (§70 are proportional to 1/3 and 1/2 power of this ratio. 



14.2. Excitation of Oblique Alfven Waves in a Kinetic Regime 

14-2.1. Cherenkov Resonance 



Using (132), (137), (]52), ( |138| ) and (132), we find a growth rate for the Cherenkov excitation of Alfven 
wave in a cold plasma: 
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with the resonant ui and k given in Table ^ for cold plasma and Table ^| for the hot plasma. In a hot plasma 
the growth rate is decreased by a factor T£. 

This growth rate is very small. Alfven waves in the limit to <C w p are almost transverse and are 
not excited effectively by the Cherenkov resonance. A strong dependence on ui and 9 corresponds to the 
increasing potential part of Alfven waves for larger uj and 9. 



14.2.2. Cyclotron Resonance 



Using (132), (137), ( ]52 ) and ( |132[ ), the growth rate for the cyclotron excitation of Alfven waves is 

r = 1 ( i49) 



4 oj ics Aj 

with the resonant frequency given in Table ^ in the cold case or Table | in the warm case. 



14.3. Excitation of the Oblique Ordinary Waves in a Kinetic Regime 

14--3.1. Cherenkov Excitation 

The Cherenkov excitation of the O mode strongly depends on the parameter \ih and the angle of 
propagation. Excitation is possible only for /j/j > 1. For fi^ > 1 the resonance occurs at w > u)q . Then 



using the polarization vector Eq. (|3^), the resonance frequency (Tables || and |J) we find from (132) 



7^ sin 2 9 



2 hr> 



A7 2 



(150) 



For the Cherenkov resonance occurs approximately at the cross-over point lOq . Using the 



polarization vector (^) , the resonance frequency @ and Eq. @) we find from Eq. (|132[) 



r = < 



7T "fc 76 

^22 " P At 3 



, sin 2 6 cos 2 6 A7 2 



(151) 



> 



Equations ( 150 ) and (151) imply that the Cherenkov excitation of the O mode is effective only if 

2uj 2 

fih ^> 1 and in the narrow angle 9 <ti —f-. This condition may be satisfied only in the outer regions of the 
pulsar magnetosphere. The growth rate of the Cherenkov excitation of the O mode in the kinetic regime is 
proportional to the density of the resonant particles. In the outer parts of pulsar magnetosphere, the density 
has decreased considerably which prevents the development of the Cherenkov instability. Numerically, it 
turns out that in the pulsar magnetosphere the kinetic instabilities may be stronger than hydrodynamic. 
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14-3.2. Cyclotron Excitation of the Ordinary Mode 



Using (136), (B2h and (132) the growth rate for the cyclotron excitation of the O wave is 



2 

Y = - ^P-re-s / 152 \ 
4 tj les cos 2 6Aj 

with the resonant frequency given in Table ^| in the cold case or Table |^ in the hot case. Here uj p . res is the 
plasma frequency of the resonant particles. The angle of emission is limited by 9 < uj p /ujb. The maximum 
growth rate, which is attained with parallel propagation, is estimated below. 



14.4. Excitation of the X Mode 



Using (135), ( p2| ) and (132) the growth rate for the cyclotron excitation of the X wave is 

r = J^%. (153) 

Using the resonant frequency (Table ||) we find 

p _ 7T ujp res ujp~fb T p _ 7rA rcs A7 b T p fl 2 (154) 
4 w|A7 A77 P lu b 

The conditions of the fast growth are 

J > 1, if (-^-] > 300 (155) 



i? c £6T 7rA rcs A 3 / 2 i? c ft / n x 3/2 



C7 2 A7 c \w B 



(156) 



Since cyclotron instability develops in the outer regions of pulsar magnetosphere, condition (156) can 
be satisfied for the regions close to the the magnetic axis with R c sa 10 10 cm. The lower streaming Lorentz 
factors increase the cyclotron instability growth rate. 

We conclude this section by the table of kinetic growth rates (Table 0) . 



15. Hydrodynamic Versus Kinetic Instabilities 

Having calculated the growth rates for the hydrodynamic and kinetic regimes of the Cherenkov and 
cyclotron instabilities, we can check whether the conditions of the corresponding regimes are satisfied. 



15.1. Cherenkov Resonance 

The condition of the hydrodynamic regime for the Cherenkov excitation is given by (^34j) with v = 
(the condition for the kinetic regime is reversed). We can distinguish two separate cases: when the scatter 
in velocity of the resonant particles is due to the the scatter in parallel velocity or to the scatter in pitch 



Resonances 


Extraordinary wave 


Ordinary wave 


Alfven wave 


Cherenkov 





r = i 

r = < 


w 2 7^ sin 2 6 ^ 1 
fc„ s c A 7 2 ' Hh > t 

• 

V22 io P A 7 2 "5" 

Hh > i 

vr rf 3 ll 2 U 2 
2 sin 2 6 cos 2 6) A 7 2 u ^ 


^ — 8 kccos6w* tan ^A 7 2 


Cyclotron 


4 w|A 7 


/ , 2 

p 7T w p,res 

4 CJres COS 2 f?A 7 


, , 2 

p TT w 6 

4 W res A 7 
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angles. In the former case condition (64) with the parallel growth rate ( |122| ) gives the following requirement 
for the hydrodynamic-type Cherenkov instability: 

7 * > 1 (157) 



which is well satisfied for the chosen plasma parameters. 

In the case when the scatter in pitch angles dominates over the scatter in parallel velocity the condition 
for the hydrodynamic type Cherenkov instability reads 

^ « (158) 

This is not satisfied. This implies that if the primary beam does not acquire any significant transverse 
gyrational energy as it propagates out in the pulsar magnetosphere, then the Cherenkov-type instabilities 
occur in the hydrodynamic regime. 

We can also verify that the condition for the kinetic growth of the beam without any scatter in pitch 
angles is not satisfied. The inverse of the condition ( |64| ) with the parallel growth rate in the kinetic regime 
(151) give the following condition for the validity of the kinetic approximation 

7fc5 < 1 (159) 



AT p 3 A 7 3 

which is not satisfied for the chosen plasma parameters. 

We conclude that the Cherenkov instability for the parallel propagation is in the hydrodynamic regime 



15.2. Cyclotron Resonance 

For the cyclotron resonances the left-hand side of ( |64| ) is dominated by the last term. For the cyclotron 
excitation of the X mode condition ( |64|) with the growth rate (|8S|) give for the hydrodynamic type instability 
to apply 

A 7 « ^%X^ = 10- ( 



(160) 



\Rns ) 

which is most probably not satisfied even in the outer regions of the pulsar magnetosphere. 



The condition for the kinetic approximation for the cyclotron excitation of the X mode follows from 
(||) and jl54D : 

a 7 » h^\K cs T p V 2 n_ (i6i) 



1p J u B 
which is well satisfied inside the pulsar magnetosphere. 

From these estimates we conclude that the cyclotron instability in the pulsar magnetosphere occurs in 
the kinetic regime. This is different from the electrostatic Cherenkov instabilities on the primary beam, 
that occur in a hydrodynamic regime. 

This difference is very important for the theories of the pulsar radio emission. The kinetic instabilities, 
in contrast to the hydrodynamic, are not suppressed by the large relativistic factor of the resonant particles. 
Thus, the kinetic instabilities are more favorable as a possible source of the pulsar radio emission. 
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It is possible to illustrate graphically the difference between the hydrodynamic regime of the Cherenkov 
instability and the kinetic regime of the cyclotron instability. On the frequency-wave vector diagram for the 
O mode (Fig. ||), the dispersion curves of the cyclotron wave in the beam oj = kv^ cos 6 — u>b/% is almost 
parallel to the dispersion curves of the excited waves in plasma at the location of the resonance. Thus, 
a small change in the velocity of the resonant particles results in a considerable change of the resonant 
frequency. This vindicates the kinetic approximation that requires a large bandwidth of the growing waves. 
In contrast, for the very large streaming 7-factor of the primary beam (so that /x, [ih 3> 1), the Cherenkov 
resonances on the O and X modes occur approximately at the cross-over frequency in a narrow frequency 
band. 



16. Conclusion 

In this work we considered normal modes and wave excitation in the strongly magnetized electron- 
positron plasma of the pulsar magnetosphere. We found the location of resonances and calculated the 
growth rates for the Cherenkov and cyclotron excitation of the O, X and Alfven waves in two limiting 
regimes of hydrodynamic and kinetic instabilities taking into account angular dependence of the growth 
rates. The main results of the paper are 

(i) Cherenkov instabilities develop in the hydrodynamic regime while cyclotron instabilities develop in the 
kinetic regime. 

(ii) Cherenkov instability on the primary beam develops on the Alfven waves in the regions close to the 



stellar surface and on the O mode in the outer regions of the pulsar magnetosphere (|73|), (115) 



(iii) Cyclotron instability can develop on all three wave branches. On the Alfven branch, the cyclotron 



instability does not develop in a very hot plasma (117) 



(iv) The typical range of angles (in the plasma frame) with the highest growth rates are 
56 ss ujp/ujg for Cherenkov excitation of the O mode 

56 w for Cherenkov excitation of the Alfven mode 
56 w ujp/ujB for cyclotron excitation of the O and X modes 

56 w 1 for cyclotron excitation of the Alfven mode 

(v) We also note, that Cherenkov instability due to the relative drift of the plasma particles can develop 
only on the Alfven mode. 

These arguments suggest that electromagnetic cyclotron instabilities are more likely to develop in the 
pulsar magnetosphere than electrostatic. 
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A. Calclations of the Resonant Integrals 

The calculations presented in this Appendix are used in all the calculations of the kinetic growth rate 
of the instabilities. 

df(p,) sfM _ 1 3 df( Pz ) 



^^(^) = __^ 7 --^J (Al) 



J dp z (kv z - wcos6») 2 f(p z )5 - ^ 
dp z (kv z -ivcosO f f(p z ) 



kc cos I ojbPz 



S(p z -Prcs) (A2) 



For lo b ^> kc this reduces to 



*^*-**h\ f ^(»j^ A (A3) 



res 



Similarly we have 



dp z {u - kv z cos6)f( Pz )5 I u — —\ = (lIM. ) (A4) 



and 



lb J \ V z 

<I,,..U - - <r -■nsflr /.,, IA ( . -^) = f^IM) ;\-,) 



lb J \ V z 



B. Relativistic Maxwellian Distribution 



We seek an appropriate expression for the relativistic one dimensional distribution. The aim of 
this Appendix is to define the relevant physical quantities measured in different systems. The relation 
obtained in the Appendix are extensively used in Section |^ when considering the properties of waves in a 
relativistically hot plasma. 



Relativistic covariant dispersion relations for plasma waves have been considered by Godfrey et al. 1974 



Melrose 1982] and others (see reference in Melrose 1982). The general expression for the frame-invariant 



distribution function is 

/(P, r) inv = ^3 exp {/i(r) - farfUu} (Bl) 

here p is a chemical potential, [3t = 1/T), T is invariant temperature, p is the momentum of the particle, 
p v is a four-momentum of the particle and U v is four velocity of the reference frame (speed of light and 
particle mass are set to unity in this Appendix). 

Next we define a flux four- vector: 

N " = f ^V/(p) inv = {n(r, t), j(r, t)} (B2) 
J 1 
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An invariant density, measured in a particular frame with the four velocity U v is then 

n = N V U„ = J ^{p u U u )f{p) im (B3) 
In particular, the invariant density in the rest frame ( with U° = {1,0,0,0}) is n° = N V U°. We normalize 



the distribution function (Bl) to the invariant density of particles in the rest frame 



n »=lww^ { -^ } (B4) 

Then, for a one dimensional distribution /(p) mv = <5(Pj_)/(p) inv / 7r (below p is a component of momentum 
along magnetic field) 

2K\{pr) 



where we introduced new variables 7 = cosh a; and 7 P = coshy and used a relation (Gradshtcin & Ryzhik 
1980| , (3.547.4)) 



/inf 
dx cosh a;exp{— /3t cosh x} = 2K\((3t) (B6) 
-inf 



The density in the frame moving with the four velocity U v — {j p , v p 7 p } ( here j p = 1/y 1 — v p ) is 

n = N° = J rfp/(p) inv - lp n (B7) 

In this work we use the distribution function normalized to the laboratory density n 

/(p) = ^/(P), f(p) = „ w exp{-/3T^^} (B8) 



There is a natural simplification of the distribution function (Bis) in the case /3t S> 1, 7 P ^> 1 (cold 
plasma streaming with large Lorentz factor). In this case the distribution is strongly peaked at 7 = 7 P so 
we can expand the distribution function, keeping terms up to the second order in 7 — 7 P : 

= nexp{-M r_M 7 -7pn = (- (7 ~ 7 ^1 (B9) 

nw mohhp I 2( 7p 2 -i) / v^A 7 I 2A 7 2 I v ; 

where we introduced A7 = VTj p and used the fact that j p 3> 1. 

In Table § we give the estimates of the moments of the relativistic Maxwellian distribution. < ... > 
implies J dp...f(p)/n, where n is the noninvariant density in the laboratory frame. The arguments of the 
Bessel functions in Table || are 1/T p . When calculating moments we used a relation 

J dxexp(- p T cosh(x))cosh(x) n = (-1)™ (BIO) 

and the asymptotic relations for the modified Bessel functions: 



K »( x ) = VT^ e ~ X ( 1 + i TT i J ( BU ) 
K (x) = - ln(z), K v {x) = \ Y(v) (f )~" x^O (B12) 



-47- 



C. Cutoff and Cross-Over Points for Parallel Propagation 

Using streaming Maxwellian distribution it is not possible to find the exact expressions for the two 
important frequencies: cutoff frequency (a limit k — > of the plasma wave dispersion) and the cross-over 
frequency (when the O mode has a vacuum dispersion relation). 



The cutoff frequency is 
and the cross-over frequency is 



2 

w cutoff 



4-7TC 2 



Aire 2 



dp f(p) 



cross— over 



7 3 {l-vf 



(CI) 



(C2) 



For the case of relativistic Maxwellian distribution (Eq. |l]) the expression for the cutoff frequency 
may be re written as 



dx 



■ exp{— (3t cosh x} 



(C3) 



cutoff 2K 1 (f3 T h P J msh 2 (x + y) 
where j p = coshy. The corresponding integrations in the case Pt S> f (cold plasma) may be performed 



using the steepest decent method. For Pt ^ 7 P the "sharply" peaked function under integral sign in ( |C3| ) 
is exp{— Pt cosh a;} so we can use expand around the point x — to obtain 

,2 



^cutoff 



7r exp{-/3 T } 
2f3 T 2rf2K 1 (f3 T ) 



if Pt > lp 



For Pt <SC 1 (hot plasma) we can make the following approximation to the exponential function: 



cxp{— Pt cosh a;} 



f if In p T < x < - In p T 
otherwise 



(C4) 



(C5) 



We find then 



^cutoff 



In 2/5 



dx 



2K 1 (p T h P y_i n2 /3 cosh 2 (o; + y) 



Tip 



for Pt~1 v > 1 
for Pt~1 p < 1 



where we used 



In 2/3 



dx 



In 2/3 cosh 2 (x + y) 



for p T < 1 



(C6) 



(C7) 



An interesting consequence of Eqs. {C4) and (C7) is that in the case of relativistically hot plasma 
streaming with very large Lorentz factor, so that 7 P > T » 1 thermal motion increases the cutoff 
frequency, while for the lower streaming Lorentz factors thermal motion decreases the the cutoff frequency. 



The calculations of the cross-over frequency (|C2|) may be done exactly 



cross — over 



2if 1 (/3 r ) 7p 



dx exp{-/? T cosh(ir - y)} _ 7 P (1 + v p ) 2 K 2 {PtW v 



cosh x (1-tanhx) 5 



2Kx(Pt) 



4 lp u; 2 p T T » f 



27 P w 2 



T< f 



(C8) 




Fig. 1 



, — Distribution function for a one-dimensional electron-positron plasma of pulsar magnetosphere. 
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Fig. 2.- 

Dispersion curves for the waves in a cold electron-positron plasma in the plasma frame for oblique propagation 
(6 = 0.5). There are three modes: Ordinary (O), Extraordinary (X) and Alfvcn. For graphic purposes the 
gyrofrequency was chosen to be = buo p . In the high frequency regime w>wg there are two subluminous 
waves with the dispersion relation oj 2 ss k 2 c 2 + 2uj 2 for the X and O modes. Both X and O modes have 
resonances at uj = lub and Alfven has a resonance at uj — V2 cos 9lo p . O mode has a cutoff at u) = ^/2oj p . O 
mode crosses the vacuum dispersion relation at uj 2 = 2oj 2 + co% sin 9 2 . 
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kc/cj p 



Fig. 3. — Dispersion curves for the waves in a cold electron-positron plasma in the plasma frame in the 
limit uo p <C uib- There are three modes represented by the dashed (O mode), solid (X mode) and long 
dashed (Alfven mode). The dotted line represents the vacuum dispersion relation. For the exact parallel 
propagation, the dispersion curves for the O mode and Alfven mode intersect. The insert in the upper left 
corner shows the region near the cross-over point ojq. 
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n=0 




W=l n=oo 



B->0 

n> 1 




Fig. 4. — CMA diagram for the X mode. The vacuum case corresponds to Z = 0. On the axis W = 
refractive index is n = 1. Resonance occurs at W = 1 ( n = oo) and reflection occurs at W = 1 + 2Z. 
Typical X waves in the pulsar magnetosphere have Z <C 1, W <C 1 and n > 1 deep in the magnetosphere. 
Arrows indicate the adiabatic tracks for the constant density and decreasing magnetic field (B — ► 0) and 
constant magnetic field and decreasing density (n e — ► 0). 



- 52 - 




Fig. 5. — CMA diagram for the and Alfvcn modes. Vacuum case corresponds to Z = 0. On the axis 
W = refractive index is n = 1. Resonances occur at n — oo and reflections occur at = 1 + 2Z and 
W = 2Z ( n = 0). The curve = 1 + 2Z corresponds to the upper hybrid wave J 2 = u% + 2uj 2 and 
the curve W — 2Z corresponds to the plasma wave ui 2 = 2ui 2 . Typical O waves (denotes by O) in the 
pulsar magnctosphere have Z« 1, 2Z < W < 1. Typical Alfven modes (denotes by A) in the pulsar 
magnetosphere have W <C 1 and n > 1. The arrows B — > 0, n e — ► and 9 — > 7r/2 indicate correspondingly 
adiabatic tracks for constant density and decreasing magnetic field, constant magnetic field and decreasing 
density and increasing angle of propagation. 
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1 2 3 4 5 

kc/o) p 

Fig. 6. — Dispersion curves for the waves in a hot electron-positron plasma in the plasma frame in the limit 
u <C u>b- Only Alfven (dotted) and O (dashed) modes are shown. The dispersion curve for the X mode is 
very similar to the cold case. For the illustrative purposes we have chosen T p = 2. The dispersion curves of 
the O and Alfven modes intersect only for parallel propagation. 
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